Equilibrium basal-plane magnetization of superconductive YNi 2 B 2 C - the influence of 

non-local electrodynamics 



O 

o 

(N 

m 

(N 
O 

o 

Oh' 
c3 



^3 

G 

O 

o 



(N 
O 



c3 



i 

c 

o 
o 



X 



J.R. Thompson 1 ^ 2 , A.V. Silhanek 3 , L. Civale 3 , K.J. Song 2 , C.V. Tomy 4 , and D.McK. Paul 5 
1 Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6061. 
2 Department of Physics, University of Tennessee, Knoxville, Tennessee 37996-1200. 
3 Comision Nacional de Energia Atomica-Centro Atomico Bariloche and Instituto Balseiro, 84-00 Bariloche, Argentina. 

4 Department of Physics, 1.1. T. Powai, Mumbai, 400076, India 
5 Department of Physics, University of Warwick, Coventry, CV4 7AL, United Kingdom. 
(submitted October 9, 2000; revised 22 February 2001) 

For a single crystal of YM2B2C superconductor, the equilibrium magnetization M in the square 
basal plane has been studied experimentally as a function of temperature and magnetic held. While 
the magnetization M(H) deviates from conventional London predictions, a recent extension of Lon- 
don theory (to include effects of non-local electrodynamics) describes the experiments accurately. 
The resulting superconductive parameters are well behaved. These results are compared with cor- 
responding findings for the case with M perpendicular to the basal plane. 



I. INTRODUCTION 

Borocarbide superconductors continue to reveal new 
and interesting features.EI This family of compounds, 
i?eNi 2 B 2 C, exhibits superconductivity for Re = Lu, Y, 
Tm, Er, Ho and Dy. The last four of these order an- 
tiferromagnetically below the Neel temperature T/v that 
ranges from 1.5 K to 10 K. The relatively high supercon- 
ducting transition temperature T c and the broad vari- 
ation of the ratio T^/T c within the family make these 
materials particularly appropriate to explore the mi- 
croscopic coexisteneeppf superconductivity and localized 
magnetic moments au 

A second remarkable feature in these compounds is the 
formation of superconducting vortex lattices rfFLL] with 
symmetries other than the hexagonal oneflQTl The pres- 
ence of these non-hexagonal lattices has been attributed 
to the effects of nonlocal electrodynamics, which arise 
when the electronic mean free path I is larger than the 
BCS zero temperature superconducting coherence length 
£o- Nonlocal electrodynamics in superconductors had 
been traditionally associated with very low values of the 
Ginzburg Landau parameter, n ~ 1. Borocarbide super- 
conductors have k values in the range of 10 to 20. How- 
ever, the availability of very clean single crystals with 
large £ permits the observation of nonlocal effects in these 
intermediate k materials. 

Nonlocality influences the equilibrium magnetic re- 
sponse of the vortex lattice in various ways. Song et 
al.E3 have shown that, when the applied field H is par- 
allel to the crystallographic c-axis of a YM2B2C crys- 
tal, the reversible magnetization in the mixed state devi- 
ates from the logarithmic dependence on magnetic field, 
M oc ln(H C 2 /-££-), which is expected from the standard 
London modelO Such deviations could be quantitatively 
accounted for within the framework of a nonlocal gener- 
alization pf. London theory, as developed by Kogan and 
Gurcvich.113 



More recently, we have measurecO M of this com- 
pound for H lying within the basal plane, and we found 
that it oscillates with angular periodicity 7r/2, a behavior 
that is incompatible with the standard London model. 
Indeed, the familiar superconductive mass anisotropy 
with components m^-, which plays a major role in the 
layered high-T c materials, is a second rank tensor. Thus, 
in the square ab basal plane of this tetragonal compound 
one has m aa = mbb, which immediately implies that the 
response within the plane should be isotropic. In con- 
trast, Kogan's nonlocal scenario! 2 ! contains a fourth rank 
tensor that breaks the basal plane isotropy and provides 
for the observed four-fold anisotropy. 

In this paper we expand the analysis of the reversible 
magnetization of the vortex system in YNi2B2C by pre- 
senting detailed measurements of M(H, T) for H in the 
basal plane of the crystal. Within the framework of the 
Kogan-Gurevich modeled, we obtain explicit expressions 
for the magnetization in the ab-plane by expanding the 
free energy appropriate for this configurationE-3 to first 
order in the basal-plane anisotropy. This approximation 
provides an excellent description of the experimental re- 
sults and, furthermore, the resulting parameters for the 
superconductor are well behaved and exhibit a remark- 
able consistency with results from band structure calcu- 
lations. These aggregate findings give unimpeachable ev- 
idence for a profound impact of nonlocal electrodynamics 
on clean, intermediate K borocarbide superconductors. 



II. THEORETICAL BACKGROUND 

Standard local London anisotropic theory provides a 
simple logarithmic field dependence for the equilibrium 
magnetization. For intermediate fields H c i <C H <C 
one hasE!l (for H parallel to the fc-th principal axis) 
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Here 77 is a constant of order unity; \ is the London 
penetration depth corresponding to screening by currents 
in the j-direction, with H || k-axis; H^ 2 = $0/2^^ is 
the upper critical field in the k-direction, with £ being 
the Ginzburg-Landau coherence length at temperature 
T . Experimentally, we will make the usual approxima- 
tion that the flux density B = H + AttM can be replaced 
by the applied field H , since the magnetization M <C H 
in all cases treated here. 

In the Kogan-Gurevich non-local formulation of Lon- 
don theoryE] there is a third independent length scale, 
the nonlocality radius p, which depends on I and T and 
also reflects the material anisotropy. It has the form 
p = X\/n, where A = (XaXbXc) 1 ^ 3 and n is the appro- 
priate component of the forth rank tensor h, given by 
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Here v is the Fermi velocity and (...) indicates aver- 
ages over the Fermi surface. The function 7(T,£), that 
contains all the temperature and mean free pa fch depen- 
dencies, has been evalulated by Kogan et alOEj 

When non-locality is important, the scale H C 2 ~ &o/£, 2 
for M is replaced by another magnetic field scale Hq ~ 
4>o/p 2 . As j(T,£) slowly decreases with increasing T, 
so does p, and consequently Hq slowly increases with 
temperature, in contrast to H c2 . One consequence of 
the theory is that the quantity jH ~ 1/£q should be 
independent of temperature, which is a prediction that 
we test later. 

In a tetragonal material h has four independent com- 
ponents, ni = naaaa, n 2 = n aa bb, 113 = n cccc and 
TI4, = n aa cc- rii or H || c-axis, the resulting expression for 
M C (H,T) isB 
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where M C = <$> /32n 2 X 2 ab , H c = $ /4:TT 2 X 2 n 2 (for a 
square FLL), and ( C (T) = m - ln{H^j] 2 H c c2 + 1), with 
both tji and r\ 2 constants of order unity. 

If H lies in the ab-plaae the analysis is more com- 
plex. In a previous studyj!3 we described the basal plane 
anisotropy in M by assuming the validity of an expression 
analogous to Eq. t^) and proposing a four-fold oscillation 
in Hq. The empirical expression so obtained successfully 
captured the basic features of the in-plane magnetic re- 
sponse, but the link between the amplitude of the oscil- 
lations and the more fundamental material parameters 
was undefined. 

According to a .subsequent generalization of the for- 
malism, develop edli-3 for the case of H in the ab-plane, 
the free energy is approximately given by 



F = Mg b B 



du 



ui y/{u + Ui){u + n 4 + d{tp)/A) 



(4) 



Here Mg b = $ /327r 2 A ah A c , m = (4 7 r 2 K 2 )^ 1 and 
u 2 = 2ttui (H c2 /B). .-(.Within this approximation, the 
very small anisotropycS in H c2 can be ignored.) Equa- 
tion (Q) holds whenever u -C 1, so that terms of order 
u 2 can be neglected. For YNigBgC we have k ~ 10, thus 
i*i ~ 2.5 x 10~ 4 and u 2 ~ 1.6 x 10~ 3 H c2 /B. In our anal- 
ysis of the basal-plane magnetization, (H c2 /B) < 40 in 
all cases (see section IV) thus u 2 < 0.05 and the approx- 
imation is valid. 

The in-plane anisotropy is accounted for by d(ip), 
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where T 2 = m c /m a is the usual mass anisotropy between 
the a and c-axes, and ip is the angle between the vortices 
and the a-axis. If d = the in-plane response would 
be isotropic and integration of Eq. (Q) would result in 
a magnetization M = —dF/dB identical to Eq. (^). 
However, for d ^ Optlic in-plane magnetization cannot in 
general be reducedt 4 ] to the form (Q). 

The integrand of (Q) can be expanded in powers of the 
variable d(ip) /A(u+rii) and integrated term by term. The 
resulting series can be differentiated with respect to B to 
obtain M ab . The leading term, M^ a , is independent of d 
and represents a large contribution to the magnetization 
that is isotropic within the plane 
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where H c lb = <&o/47T 2 A 2 ri4 (Also in this case, the numer- 
ical factors are valid for a square FLL). The term C, ab 
arises from the core contribution to the free energy (not 
included in Eq. (0)); it is analogous-, to. C c in both form 
and origin, as described previouslyliirL 4 ] Eq. has the 
same functional form as Eq. (|J), but the anisotropy be- 
tween the c-axis and the plane is reflected in both the 
prefactor Mq and the field scale Hq. 

The second term in the expansion, Mf>, is linear m 
d(tp) and accounts for the in-plane fourfold anisotropy in 
M ab to leading order: 
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In the above scenario H c2 tends toward zero as T — > T c , 
while Hq increases with temperature. These differing 
temperature dependencies mean that for clean materi- 
als, the non-local expressions Eqs. || and || reduce to 
the local form, Eq. [I], as T approaches T c , while Mf b 
from Eq. [7] vanishes. Thus, the non-local theory pre- 
dicts that the equilibrium magnetization in a clean sam- 
ple should vary logarithmically with field near T c , but 
deviate progressively from logarithmic behavior at low 
temperatures. Furthermore, as materials become dirtier 
so that p becomes shorter and Hq increases, the non-local 
expressions reduce to the local form at all temperatures. 
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Our purpose is to analyze the in-plane magnetization 
using Eqs. f| and [?]. We thus need to estimate the impor- 
tance of the various terms in-the expansion. According 
to electron band calculations,!!!! for YN12B2C we have 
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Based on those values we can calculate the relations 
between all the components of n using Eq. (||). In section 
IV we will compare these band calculation estimates with 
our experimental results. In particular, we are interested 
in the dimensionless prefactor d/8114 that sets the order 
of magnitude in Eq. m. It is also useful to split d in two 
parts, d — di + d2 sin (2ip). Using Eq. we obtain 

A- « 0.23 A « -0.13 (9) 

8ri4 8ri4 

In section IV we will make use of these estimates to 
gain an idea of the goodness of our approximations. For 
comparison, we can calculate the equivalent values for the 
similar material LUN12B2.C, using the Fermi-surface av- 
erages given previously.tLa The results are d\/&n$ w 0.42 
and ^2/8714 ss —0.20. The larger values for the Lu-based 
system mean that the first-order expansion in Eqs. 6-7 
is more accurate for the yttrium-based compound, while 
the amplitude of oscillation in the basal plane magneti- 
zation can be larger for LUM2B2C. 

III. EXPERIMENTAL ASPECTS 

The YM2B2C single crystal was grown by a high tem- 
perature flux method using N12B flux, using isotopic n B 
to reduce neutron absorption in complementary scatter- 
ing studies. The 17 mg crystal is the.same as that used in. 
previous investigations by Song et al£j and Civale et al.E3 
It is a slab of thickness t ~ 0.5mm in the c-axis direction, 
with a mosaic spread of less than 0.2°, as determined by 
neutron diffraction. In the basal plane the shape is ap- 
proximately elliptical with principal axes of length ~ 2.0 
and 2.5mm, which approximately coincide with the two 
equivalent (110) axes of the tetragonal structure. 

Magnetic studies were conducted in a SQUID-based 
magnetometer Quantum Design MPMS-7 equipped with 
a compensated 70 kOe magnet. Normally, scan lengths of 
3cm were used. The crystal was glued onto a thin Si-disk 
and mounted in a Mylar tube for measurement with the 
magnetic field H applied in the basal plane, along either 
the [100] or the [110] axis, with an accuracy better than 



3° . For both orientations, the magnetization M is paral- 
lel to H by symmetry, thus only the longitudinal compo- 
nent measured by the magnetometer needs to be consid- 
ered. This is also true with H | [001], as was the case for 
the previous data that will be included in our analysis. 
Then, from now on we will ignore the vector nature of M 
and denote it simply as M. The diamagnetic moment of 
the addenda (silicon disk plus glue) , which was measured 
separately, was linear in iJ, isotropic, and non-hysteretic. 
This signal, msi = (—5.4 x 10~ 9 emu/Oe)H , was always 
small compared with the moment of the crystal, and was 
subtracted from all the data prior to any further analysis. 

In the mixed state, hysteresis loops M(H, T) were mea- 
sured. The maximum H was 65fcOe in all cases. Mea- 
surements were also conducted in the normal state at 
temperatures up to 300-fT, in order to correct for the 
normal state background moment. It is worth noting 
that the magnetization is small, compared to the applied 
field, in all cases considered here. Thus demagnetizing 
efects are negligible: for H_|l a6-plane, we obtain from 
the Meissner state responsefcj that the effective demag- 
netization factor D w 0.1. Then in the mixed state, the 
effective field H e f t = H app u e d — AirDM differs from the 
applied field by less than 1 % is the worst case and we 
need to consider only the magnetizing field "iJ." Fur- 
thermore, for comparison with theoretical expressions, 
we can approximate the flux density B = H + AttM by 
H with an accuracy of a few percent and generally better. 

The superconductive transition temperature, mea- 
sured in a small applied field, was T c = 1A.5K . Measure- 
ment of the electrical resistivity using a van der Pauw 
method gave an electrical resistivity of Afifl — cm at 20K 
and a residual resistance ratio of 10, yielding an elec- 
tronic mean free path £ « 300 A. Using this and values 
of H c2 (for H 1 1 c-axis), Song et al. deduced the values 
£0 = 120 A for the BCS coherence length at T = and 
k ~ 10. Other superconducting parameters for this com- 
pound are collected in Table I. 

IV. RESULTS AND DISCUSSION 

For magnetic field orientations H j| [100] and H | 
[110], we measured isothermal magnetization loops 
M (H) at temperatures T — 3 to 1AK, in IK intervals, 
and also at 18K, slightly above T c . Three of those loops, 
with H || [110], are shown in Fig. 1, at temperatures 
T = 5K, 12K and 18K, i.e., well below, near, and just 
above T c . The magnetic response of this material in the 
superconducting mixed phase is slightly irreversible, re- 
flecting weak pinning of flux lines. Now, as the only 
source of magnetic hysteresis is vortex pinning, M(H) 
becomes reversible for H > H C 2 (T) . Based on the Bean's 
critical state model, we calculated the equilibrium mag- 
netization as the average M eq (H) = [M T + M l ] /2 of 
the magnetization values measured in the field-increasing 
and field-decreasing branches of the loop, respectively. 



3 



The result for T — 5K is shown as a dashed line in Fig. 
1. The identification of M eq (H) with the average be- 
tween the branches of the loop increases in accuracy as 
the width of the hysteresis loop \M' — MH decreases. 
For each temperature, we disregard the low field data 
where large hysteresis introduces a significant uncertainty 
in the determination of M eq . Experimentally, as the 
field decreases, the width of the M(H) loop increases 
continuously and smoothly. Then at some temperature- 
dependent low field, M changes abruptly as it reaches a 
minimum and starts to increase. This feature is visible 
in Fig. 1 for T — 5 and 12 K, for example. Below this 
field, the hysteresis grows suddenly. We compute the av- 
erage magnetization starting at a field slightly above this 
minimum. 

It is apparent in Fig. 1 that the magnetization M eq 
does not vanish above H C 2 and above T c . This indicates 
the presence of a normal state contribution M ns (H, T) 
to the magnetization (recalling that the signal from 
the sample holder has been already subtracted). Pure 
YNi 2 B 2 C has no localized moments, thus in the nor- 
mal state it is expected to exhibit a linear and nominally 
temperature-independent paramagnetic (Pauli and Van- 
Vleck) susceptibility xo- However, close inspection of the 
data indicates that this intrinsic term cannot account for 
the entire normal state signal. Indeed, M ns (H,T) is not 
linear in H and it grows as T decreases, thus pointing to 
the presence of localized moments. The small magnitude 
of the signal suggests, on the other hand, that this contri- 
bution arises from magnetic impurities. To confirm this, 
we measured the temperature dependence of M ns (H, T) 
from 16K to iQQK at several fixed fields. The results for 
H = 1 and 5kOe are shown in the inset of Fig. 1. As 
the localized moments are very dilute, no magnetically 
ordered phase should appear and one can expect a Curie 
Law dependence. The solid lines in the inset are fits to 
M ns /H = [xo + CJT\. The Curie term corresponds to a 
rare earth impurity content of ~ 0.1 at % relative to yt- 
trium, most likely contaminants in the yttrium starting 
material. 

To isolate the magnetization associated with the vortex 
state, it is necessary to remove the paramagnetic back- 
ground. The normal state magnetization M ns is well 
described by the expression 

M ns = xvH + M sat Bj (^?) (10) 

where xo is t ne temperature independent susceptibility 
obtained as explained above and Bj is the Brillouin func- 
tion for effective angular momentum J. For small values 
of the argument H/T, the Brillouin function leads to the 
Curie susceptibility, of course. 

The magnetic signal arising from the localized mo- 
ments is clearly observed in Fig. 2, where [M eq — XoH] is 
plotted as a function of H/T. The figure includes exper- 
imental data from temperature sweeps at three values of 
H, and from isothermal loops at temperatures down to 



hK. At each T, the onset of the superconducting signal 
below H c2 (T) is clearly observed in the nearly vertical 
trace of data below the envelope curve. It is also ap- 
parent that for H > H C 2(T) all the data at different T 
overlap on a single curve arising from the paramagnetic 
contribution of the magnetic impurities. The solid line 
is a nonlinear fit to Eq. [lO] that yields parameter val- 
ues J = 3/2 and g = 6. The saturation magnetization 
M sat — 0.92G corresponds to ~ 7 x 10 -4 per formula 
unit content of rare earth impurities, most likely trace 
contaminants in the yttrium starting material. At this 
concentration, the magnetic impurity ions are isolated, 
but strongly influenced by the crystal field of the host, 
producing an anisotropic magnetic response in the nor- 
mal state. The susceptibility is smaller for £f||c-axis, 
as found for i?eNi2B*C crystals with Re = Tb, Dy, 
and Ho, but not Tm.E3 These qualitative similarities 
show that (a mixture of) any of several rare earth impu- 
rities can generate the observed paramagnetic response. 
The central and important point here, however, is that 
the empirical fit to Eq. [l^ provides a precise description 
of the normal state signal, so that the superconductive 
magnetization can be isolated for analysis. 

The resulting (background-corrected) superconducting 
state equilibrium magnetization M — M eq — M ns is 
shown in Fig. 3, as a function of magnetic field H applied 
along the [110] axis of the crystal. Results are shown for 
temperatures 3K < T < 14K in intervals of IK. Qual- 
itatively, the curves for temperatures near T c are linear, 
showing the dependence on ln(B) as predicted by tradi- 
tional local London theory. At lower temperatures, how- 
ever, the increasing curvature visibly signals a progressive 
departure from local London behavior. 

For a quantitative analysis, we fit the low temperature 
data to the non-local relation, Eq. (^), varying the pa- 
rameters M$ b , H% b , and ( ab . These fits describe the low 
temperature experiments very well, as shown by the solid 
lines in Fig. 3. The resulting values of Mfi b , Hg b , and Q ab 
as a function of T (up to T = 1QK) are shown in Figs. 
4(a), (b) and (c) respectively. At higher temperatures, 
T > 11K, the system closely approximates local London 
behavior, with M cx ln(B) and we analyze these data 
using the local expression Eq. ([[]) . Numerically, the data 
close to T c have a very shallow minimum as a function 
of Ho and it becomes meaningless to fit an (essentially) 
straight line with three parameters. The fits to the lo- 
cal London relation are shown in Fig. 3 as dashed lines, 
and the values of Mq(T) obtained by this procedure are 
also included in Fig. 4(a) (open symbols). In the non- 
local analysis, H c i was not a fitting parameter. It was 
determined independently, as explained below. 

Before we discuss the results shown in Fig. 4, we must 
analyze what is the error involved in taking M ab sa Af?* , 
disregarding terms of higher order. To estimate the con- 
tribution M£^ given by Eq. |?j, we must rely on the numer- 
ical values CJ9|) obtained from band calculations. We see 
that, if at least the signs in (g) are correct, Mf 6 minimizes 
for ip = 45°. For that reason we chose to apply Eq. ^| to 
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the [110] data and not to the [100] orientation. Numer- 
ically, we have d(ip = 45°)/8n 4 = [di + d 2 ) /8n 4 w 0.10. 
Thus, Eqs. @ and (@) indicate that Mf b < Mf s b 
over most of the field range (for instance, at T = 3K 
and B = 2kG, M? b /M$ ~ 0.02, and the same ra- 
tio is obtained for 6K and lOkG). We can also com- 
pare the logarithmic field derivatives of Af, which give 
the slopes in Fig. 3. For B <C Ho, to leading or- 
der we obtain dM^JdhxB w Mg b and dMf/dlnB w 
0.10 {2B/H ) M ab , thus for B = 2fcG the error in the 
slope produced by disregarding the second term in the 
expansion is less than 0.5%. 

The above estimates indicate that by using Eq. (^) to 
analyze the magnetization in the [110] direction we are 
introducing an error in the determination of M^ h and 
i?Q b of at most a few percent, which is less than the 
experimental noise observed in Figs. 4(a) and (b). So, 
it is meaningful to discuss the results shown in Fig. A. 
For comparison, the corresponding results for the casecHI 
with H || [001] are included, too. 

Figure 4(a) shows that M varies linearly with T 
near T c and extrapolates to zero at T c , consistent with 
Ginzburg-Landau theory. From Eq. (|l|), the ratio 
M§/Mg b = \ c /\ab = \/{rn c /m ab ) = T. As expected, we 
find that this ratio is rather independent of temperature, 
r = 1.13 ± 0.02. This result is very comparable with 
the experimental value of 1.16 for LUN12B2C obtained 
by Metlushko et al.E^I This experimental value can also 
be compared with the band structure calculation, since 
T = \/(( v a) /( v c))- The resulting value, 1.01, is much 
smaller than that deduced experimentally. Returning to 
the experimentally determined parameters, we can also 
calculate A C (T), A afc (T) and A(T). The results at T = 3K 
are shown in Table 1. 

The next frame, Fig. 4(b), shows Ho(T). Qualita- 
tively, Ho for both orientations is constant at low tem- 
perature, then increases as T increases., consistent with 
the theoretically predicted behavior£3ll3 In section II 
above, we noted that the non-local theory predicts that 
the quantity i?o 7 should be independent of temperature. 
PreviouslytLEl we found that £0 = 12 nm; this value agrees 
within experimental error with that calculated from the 
BCS expression £ = h{v f) / {k 2 e.^ 1 k B T c ) = 11.6 bhi, us- 
ing (vf) = 3.87 x 10 7 cm/sec from band structure, til with 
e 7 w 1.78. Also, from the electrical resistivity, we have 
that £ w 30 nm for this crystal and therefore evaluate 
the impurity parameter j(T) with £o/£ = 0.3. The re- 
sults for i?o7 are shown as open symbols in Fig. 4b, for 
the two orientations of magnetic field. The constancy 
of the product provides further solid evidence that non- 
local electrodynamics strongly modify the superconduc- 
tive properties of clean borocarbides. From the Ho data 
we can extract the nonlocality radius in both orienta- 
tions, p ab (T) = ($ /47T 2 iJ§) 1/2 (see Song et al.0), and 
analogous expression for p c (T). We can also calculate 
two of the four independent components of the tensor n, 
namely n 2 — (pab/ty 2 and m = (p c /X) 2 . All the numer- 



ical values for T = 3K are listed in Table 1. 

Figure 4(c) shows the temperature dependence of the 
fitting parameter £. Qualitatively, the data reflects the 
fact that C ~ —^(Ho/Hcs). As T increases, Ho/H C 2 
varies little at low temperatures, but then becomes 
larger, due to the differing dependencies of Ho{T) and 
H c2 (T), thus C(T) decreases. 

The last frame, Fig. 4d, shows H c2 (T) in both ori- 
entations, obtained (a) by extrapolating the isothermal 
magnetization in the superconductive state to M = 
and (b) by locating the field at which the magnetic hys- 
teresis disappears. These results were introduced as fixed 
parameters in Eq. (^), as already mentioned. It is ev- 
ident that H C 2(T) curves upward near T c . This differs 
from simple Ginzburg-Landau behavior, bt)fjf consistent 
with other observations on this material. Also, little 
anisotropy is evident between the c-axis and the basal 
plane values. The mass anisotropy obtained from the Mo 
data (Fig. 4(a)) would imply an experimentally resolv- 
able, difference in H c2 . Previously, Johnson-Halperin et 
al.t£l also found a nearly isotropic response in single crys- 
tal YNI2B2C. The reason for the discrepancy between 
the anisotropy derived from Mo and H c2 remains un- 
clear. Interestingly, LuNi2B2C, an isostructural, non- 
magnetic borocarbide superconductor with similar T c , 
exhibits greater anisotropy H c2 , both from the c-axis 
and within the basal plane .113 

Next we consider the anisotropy of the magnetization 
within the basal plane. We have previously shownt3 that 
this quantity exhibits a four-fold periodicity, which can- 
not be accounted for within the local London theory. 
In Fig. 5 we show the amplitude of the oscillation in 
basal plane magnetization, SM = M[110] — A/ [100], as a 
function of field H. The experimental data for T = 7K 
are reproduced from Civale et al.Ej, where details of the 
experimental procedure and data analysis are provided. 
Here we analyze these data using Eq. (0), where only 
the second term, proportional to d 2 , contributes to the 
oscillatory behavior. In this expression, we have already 
values for H c2 and -ffp b , so there is a single unknown, 
d 2 /8ii4, available for fitting the data. The result is shown 
as a solid line in Fig. 5. The value of the sole fitting pa- 
rameter is ^2/877.4 = —0.136, which we discuss below. 

Let us now summarize the experimentally determined 
parameters related to the material anisotropy. Theoreti- 
cally, five independent parameters {ni,n 2 , 113,114 and T) 
are necessary to describe fully the angularly dependen- 
cies. Experimentally we have obtained three of these 
quantities, n 2 , 114 and T, plus the combination d 2 /8n4 
that allows us to obtain the value for n\ using Eq. 5. 
Thus, we have obtained four out of the five indepen- 
dent parameters, as shown in Table I. For comparison 
with band structure calculations and to eliminate poorly 
known prefactors, it is convenient to consider ratios rel- 
ative to ri4, as shown in Table II. The experimental and 
calculated ratios agree within ~ 15 %, which is quite 
reasonable given the approximations. Lastly we consider 
the sole fitting parameter d 2 /8ri4 ~ —0.14 that deter- 
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mines the amplitude of the oscillations of M in the basal 
plane. This quantity has the correct sign and it lies in re- 
markable numerical agreement with the band structure 
value, -0.13. The excellent description of the data, us- 
ing parameter values measured independently here (M , 
H c2 , and H ) and with ^2/8/14 almost coinciding with 
the calculated value, gives considerable confidence in the 
fundamental correctness of the non-local description. 



V. CONCLUSIONS 

We have measured the equilibrium magnetization in 
the basal plane of clean, superconducting YNi 2 B 2 C, tak- 
ing careful account of background effects. Influences of 
non-local electrodynamics are very evident well below T c , 
with the effects becoming washed out at higher tempera- 
tures as expected. The magnetic-field-dependence of the 
magnetization and its oscillation amplitude are well de- 
scribed by a generalization of London theory to include 
non-local influences. The results of this analysis have 
been compared with a corresponding study of the c-axis 
magnetization, thereby providing the value T = 1.13 for 
the mass anisotropy of this material. Furthermore, the 
experimental values for the material parameters agree 
overall very well with those deduced from the band struc- 
ture. In summary, these clean, non-magnetic borocarbidc 
superconductors display a rich variety of physical phe- 
nomena, requiring both 2-nd and especially 4-th rank 
tensors to describe their macroscopic vortex state prop- 
erties. 



VI. ACKNOWLEDGEMENTS 

We wish to acknowledge useful discussions with V.G. 
Kogan, M. Ycthiraj and D. K. Christen and to thank 
W.E. Pickett for his unpublished band structure results. 
A portion of the work of JRT was supported by the Sci- 
ence Alliance at The University of Tennessee, Knoxvillc. 
Research at ORNL was sponsored by the U.S. Depart- 
ment of Energy under contract DE-AC05-00OR22725 
with the Oak Ridge National Laboratory, managed by 
UT-Battelle, LLC. Research at CAB-IB was partially 
supported by ANPCyT, Argentina, PICT 97 No. 01120, 
and by CONICET, Argentina, PIP No. 4207. A.V.S. 
would like to thank the CONICET for financial support. 



VII. REFERENCES 



1 P.C. Canfield, P.L. Gammel and D.J. Bishop, Physics To- 
day 51, 40 (1998). 

2 U. Yaron, P.L. Gammel, A. P. Ramirez, D.A. Huse, D.J. 
Bishop, A.I. Goldman, C. Stassis, P.C. Canfield, K. 
Mortensen, and M.R. Eskildsen, Nature 382, 236 (1996). 

3 H. Eisaki, H. Takagi, R. J. Cava, B. Batlogg, J. J. Kra- 
jewski, W. F. Peck, Jr., K. Mizuhashi, J. O. Lee, and S. 
Uchida, Phys. Rev. B 50, 647 (1994). 

4 M. R. Eskildsen; K. Harada; P. L. Gammel; A. B. Abra- 
hamsen; N. H. Andersen; G. Ernst; A. P. Ramirez; D. J. 
Bishop, K. Mortensen, D. G. Naugle, K. D. D. Rathnayaka, 
and P. C. Canfield, Nature 393, 242 (1998). 

5 P. L. Gammel, B. P. Barber, A. P. Ramirez, C. M. Varma, 
D. J. Bishop, P. C. Canfield V. G. Kogan, M. R. Eskildsen, 
N. H. Andersen, K. Mortensen, and K. Harada, Phys. Rev. 
Lett. 82, 1756 (1999). 

6 K. N0rgaard, M. R. Eskildsen, and N. H. Andersen, J. 
Jensen, P. Hedegard, and S. N. Klausen, and P. C. Canfield, 
Phys. Rev. Lett. 84, 4982 (2000). 

7 M. Yethiraj, D.McK. Paul, C.V. Tomy and E.M. Forgan, 
Phys. Rev. Lett. 78, 4849 (1997). 

8 Y. De Wilde, M. Iavarone, U. Welp, V. Metlushko, A.E. 
Koshclev, I. Aranson, GW. Crabtree and P.C. Canfield, 
Phys. Rev. Lett. 78, 4273 (1997). 

9 D.McK. Paul, C.V. Tomy, CM. Aegerter, R. Cubitt, S.H. 
Lloyd, E.M. Forgan, S.L. Lee and M. Yethiraj, Phys. Rev. 
Lett. 80, 1517 (1998). 

10 K.J. Song, J.R. Thompson, MYethiraj, D.K. Christen, 
C.V. Tomy and D. McK. Paul, Phys. Rev. B 59, R6620 
(1999). 

11 V.G. Kogan, M.M. Fang, and S. Mitra, Phys. Rev. B 38, 
11958 (1988). 

12 V.G. Kogan, A. Gurevich, J.H. Cho, D.C. Johnston, Ming 
Xu, J.R. Thompson and A. Martynovich, Phys. Rev. B 54, 
12386 (1996). 

13 L. Civale, A.V. Silhanek, J.R. Thompson, K.J. Song, C.V. 
Tomy and D.McK. Paul, Phys. Rev. Lett. 83, 3920 (1999). 

14 V.G. Kogan, S.L. Bud'ko, P.C. Canfield and P. Miranovic, 
Phys. Rev. B 60, R12577 (1999). 

15 V. G. Kogan, P. Miranovic, and D. McK. Paul, in The Su- 
perconducting State in Magnetic Fields: Special Topics and 
Trends, edited by C. A. R. Sa de Melo, Directions in Con- 
densed Matter Physics, vol 13 (World Scientific, Singapore, 
1998). 

16 E. Johnston-Halperin, J. Fiedler, D. E. Farrell, Ming Xu, 
B. K. Cho, P. C. Canfield, D. K. Finnemore, and D. C. 
Johnston, Phys. Rev. B 51, 12852 (1995). 

17 W.E. Pickett (unpublished). 

18 V.G. Kogan, M. Bullock, B. Harmon, P. Miranovic, Lj. Do- 
brosavljevic-Grujic, P.L. Gammel, and D.J. Bishop, Phys. 
Rev. B 55, R8693 (1997). 

19 V. Metlushko, U. Welp, A. Koshelev, I. Aranson, G. W. 
Crabtree, and P. C. Canfield, Phys. Rev. Lett. 79, 1738 
(1997). 

20 B. K. Cho, Physica C 298, 305 (1998). 

21 S. V. Shulga, S. -L. Drechsler, G. Fuchs, K. H. Mueller, K. 
Winzer, M. Heinecke, and K. Krug, Phys. Rev. Lett. 80, 
1730 (1998). 

22 G. Hilscher and H. Michor in Studies of High Temperature 
Superconductors, vol. 28, edited by A. V. Narlikar (Nova, 



G 



New York, 1999), pp. 28-74 and references therein. 



TABLE I. Superconducting parameters 



A = 9504 
X ab = 9901 
p ab = 31.4A 
ni = 6.95.10" 

n 3 = - 



r= 1.13 ±0.02 

A c = 880A 
p c = 34.4A 
n 2 = 1.04.1GT 3 
n 4 = 1.25.1GT 3 



TABLE II. Average Fermi velocities 




exp. 


band calc. 


n\/n& 


5.56 


4.78 


112/114 


0.83 


0.73 


n 3 /n 4 




3.2 



Figure captions 

Fig. 1. The magnetization of single crystal YM2B2C 
at three temperatures, with magnetic field in the basal 
plane, H || [110] axis. The figure illustrates the lim- 
ited irreversibility in the superconductive state below 
T c = 1A.5K, and a paramagnetic background in the nor- 
mal state. Inset: Curie analysis of the normal state sus- 
ceptibility in the temperature range T = 18 — 300A". 

Fig. 2. Total equilibrium magnetization M eq minus 
XoH (see text) versus H/T obtained from isothermal field 
sweeps at several T and from temperature sweeps at three 
values of H, as indicated in the figure. The solid line is 
a Brillouin function (Eq. [Io|) fitted to the paramagnetic 
impurity-background data. The nearly vertical traces of 
data show the entry into the superconductive state. 

Fig. 3. The equilibrium magnetization M eq in the su- 
perconductive state with H || [110] axis, plotted versus H 
(logarithmic axis). Discrete symbols are data measured 
at temperatures of 14, 13, 12, 3 K; straight (dashed) 
lines show conventional, local London behavior near T c . 
A nonlocal generalization of London theory (Eq. ^, solid 
lines) accounts well for the pronounced deviations from 
local behavior at lower temperatures. 

Fig. 4. Superconductive parameters of YM2B2C with 
magnetic field H || [001] or H 1 1 [110]. (a) The magnitude 
of the equilibrium magnetization Mo, where solid sym- 
bols denote results from the non-local analysis and open 
symbols come from a conventional London analysis, (b) 
The field scale H and the (theoretically constant) quan- 
tity H j (where 7 = impurity parameter; see text), (c) 
The fitting parameter (. (d) Estimates of the upper crit- 
ical field H C 2- 

Fig. 5. Amplitude of the oscillation in basal plane 
magnetization, 6M = M[110] — M[100], plotted versus 
H on a logarithmic axis. Solid line shows a fit to the 
oscillatory term in Eq. , with Hg b = 48 kOe and 
= 33 kOe; see text. 
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d 2 /8n 4 = 


0.136(2) 


M ab = 


5.29 G 


ij ab _ 

H o - 


48kOe 


tt ab _ 


33kOe 


R 2 = 


0.9949 
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